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Abstract 

In light of recent lattice results for the light quark masses, m s and m u + m^, 
we re-examine the use of sum rules in the extraction of these quantities, and 
discuss a number of potential problems with existing analyses. The most impor- 
tant issue is that of the overall normalization of the hadronic spectral functions 
relevant to the sum rule analyses. We explain why previous treatments, which 
fix this normalization by assuming complete resonance dominance of the contin- 
uum threshold region can potentially overestimate the resonance contributions 
to spectral integrals by factors as large as ~ 5. We propose an alternate method 
of normalization based on an understanding of the role of resonances in Chiral 
Perturbation Theory which avoids this problem. The second important uncer- 
tainty we consider relates to the physical content of the assumed location, s , of 
the onset of duality with perturbative QCD. We find that the extracted quark 
masses depend very sensitively on this parameter. We show that the assumption 
of duality imposes very severe constraints on the shape of the relevant spectral 
function in the dual region, and present rigorous lower bounds for m u + as a 
function of So based on a combination of these constraints and the requirement of 
positivity of ps{s). In the extractions of m s we find that the conventional choice 
of the value of So is not physical. For a more reasonable choice of Sq we are 
not able to find a solution that is stable with respect to variations of the Borel 
transform parameter. This problem can, unfortunately, be overcome only if the 



hadronic spectral function is determined up to significantly larger values of s 
than is currently possible. Finally, we also estimate the error associated with the 
convergence of perturbative QCD expressions used in the sum rule analyses. Our 
conclusion is that, taking all of these issues into account, the resulting sum rule 
estimates for both m u + and m s could easily have uncertainties as large as a 
factor of 2, which would make them compatible with the low estimates obtained 
from lattice QCD. 
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1 INTRODUCTION 



The recent lattice results for the light quark masses [1], m u + m d = 6.8 ± 0.8 ± 0.6 MeV 
and m s = 100 ± 21 ± 10 MeV in the quenched approximation, and the even smaller values 
m u + m d = 5.4 ± 0.6 ± 0.6 MeV and m s = 68 ± 12 ± 7 MeV, for the n f = 2 flavor theory, (all 
evaluated in the MS scheme at fi = 2 GeV), appear to be significantly smaller than results 
obtained from sum rule analyses. The most recent and complete sum rules analyses are 
those of (i) Bijnens, Prades and deRafael (BPR), which yields m u + m d (p = 2 GeV) = 9.4± 
1.76 MeV [3], and (ii) Chetyrkin, Pirjol, and Schilcher (CPS) which gives m s = 143 ±14 MeV 
[4]. We have translated the original values, m u + m d = 12±2.5 MeV and m s = 203±20 MeV, 
quoted at a — 1 GeV, to a = 2 GeV using the renormalization group running and the 
preferred value Aq CD = 300, 380 MeV used respectively in the two calculations. (The 
analysis by CPS is an update of that by Jamin and Miinz (JM) [5], however, since the 
approach and techniques are the same, we will refer to their work jointly by the abbreviation 
JM/CPS.) The sum rule results, thus, lie roughly 1 — 2a above the quenched results. The 
difference between the sum rule and the n/ = 2 lattice estimates, however, is large and, we 
feel, significant enough to warrant scrutiny. Both the lattice and sum rules approaches have 
their share of systematic errors. A recent review of the lattice results is given in [2]. Here 
we present a re-evaluation of the sum rules analyses. 

The issues in the sum-rule analyses that we shall concentrate on are the convergence of 
perturbative QCD (pQCD) expressions, the choice of s - the scale beyond which quark- 
hadron duality is assumed to be valid, and the normalization of resonance contributions in 
the ansatz for the hadronic spectral function for s < s . 

The first issue is important because both the a s and a 2 s corrections to the 2-point corre- 
lation functions used in sum rule analyses are large. This issue has been analyzed in detail 
by CPS for the extraction of m s ; therefore, we shall only comment on it briefly for the case 
of m u + m d . 

The second point is important because, as we will show below, it turns out that the 
extraction of the quark masses, in particular that of m u + m d , is very sensitive to the choice 
of so • This is illustrated by deriving lower bounds on m u + m d associated with the positivity 
of P5(t), and by investigating trial spectral functions. Ideally one would like to pick s large 
enough so that pQCD, to the order considered, can be shown to be reliable. Unfortunately, 
for larger Sq, the hadronic spectral function receives contributions from an increasing number 
of intermediate states, and hence becomes increasingly hard to model. We discuss the 
uncertainties introduced by a compromise choice of sq. In the extraction of m s by JM/CPS 
we argue that an artificially large value of Sq has been used. For a more reasonable value 
of s we are n °t able to find an estimate for m s that is stable under variations of the Borel 
transform scale u. 

The third issue arises because the continuum part of the hadronic spectral function is 
typically represented as a sum-of-resonances-modulation of a continuum form, the overall 
normalization of which is fixed by assuming complete resonance dominance of the spectral 
function near continuum threshold. This turns out to be potentially the most important 
issue. We, in fact, show in the case of the vector 2-point function, for which experimental 
information on the spectral function is available in the resonance region, that an analogous 
extrapolation from threshold to the p meson peak would lead to an overestimate of the 
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spectral function in the resonance region by a factor of ~ 5. We then explain the origin 
of this problem from the point of view of the existing phenomenological understanding of 
how resonance contributions enter the expressions for low-energy observables as computed in 
Chiral Perturbation Theory (xPT). Based on this understanding, we propose an alternate 
method for normalizing the spectral function in the resonance region which requires as input 
only the expression obtained from %PT to 1-loop order, in the near-threshold region. We then 
employ this method in a re-analysis of the only sum rule treatment for which the relevant 
%PT expression is known, namely that of the correlator of the product of two divergences of 
the strangeness-changing vector current (as used by JM/CPS to obtain the estimate quoted 
above for m s ) and show that the traditional method of normalization leads to a significant 
overestimate of m s . 

We find that the size of the corrections suggested by our consideration of the above issues 
can easily lower the sum rule estimates for both m u + m d and m s by a factor at least as large 
as two. In particular, using the corrected normalization for the hadronic spectral function in 
the JM/CPS analysis alone would lower the extracted value of m s by almost exactly a factor 
of 2. Such a change would make the lattice and sum rule estimates consistent. Lowering 
both estimates by roughly the same factor would, moreover, preserve agreement of the ratio, 
r = 2m s /(m u + m d ), with that predicted by x?T. 

The paper is organized as follows. In order to make it self-sufficient, and to introduce 
notation, we reproduce the necessary details from Refs. [3] and [5] in Sections 2 and 7. The 
convergence of pQCD is discussed in Section 3. In Section 4 we derive lower bounds on 
m u + ma-, as a function of sq, using the positivity of the relevant spectral function, p 5 . In 
Section 5 we illustrate the potential sensitivity of the extracted value of m u + m d to the 
choice of sq by considering a number of plausible trial spectral functions. The important 
issue of the overall normalization of the hadronic spectral function is investigated in Section 6 
using the vector current case as an illustrative example. Based on the lessons learned from 
the vector channel, a re-analysis of the JM/CPS estimate of m s is presented in section 7. 
Finally we end with some conclusions in Section 8. 

2 FINITE ENERGY SUM RULES 

The standard starting point for the extraction of the light quark mass combination, m u + m d , 
is the Ward identity relating the divergence of the axial current to the pseudoscalar density, 

d*A&(x) = (m d + m^x^h^hq^) (i) 

where q = (u, d,s), and the projections ± = (A x ± iA 2 )/2 pick out states with the quantum 
numbers of the ti ± . This relation implies, for the two-point function of the product of two 
such divergences, that 

^ 5 (g 2 ) = ij rf 4 xe^(0|r{9M( t -)(x),9M( / +)(0)}|0), 

(m d + m u ) 2 ij d 4 xe iq - x (0\T{P^(x),P( + \0)}\0} . (2) 

The idea of the standard analysis [3, 6, 7, 8] is then to consider the finite energy sum rules 
(FESR) generated by integrating products of the form t n ^^{t) over the contour shown in 
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Fig. 1. For n negative the result involves ^5 or its derivatives at t — 0, while for n greater 
than or equal to zero, the result is zero. For sufficiently large radii, s, of the circular portion 
of the contour, the pseudoscalar two-point function, and hence also its line integral over the 
circle, can be evaluated using perturbative QCD. Taking the resulting expressions to the 
RHS's, one obtains FESR's for the moments of the spectral function p 5 (t) = ^lm^ 5 (t), on 
the interval (0, s), for example[3], 

[dtp,{t) = ^K( S )+m,( S )] 2 ^ |l + J R 1 ( s ) + 2^i^|; (3) 

j o S dttp 5 (t) = ^K( S )+m,( S )] 2 ^|l + i? 2 ( S )-^^^}, (4) 

where m(s) is the running mass evaluated at the scale s, Ri(s) and #2(5) contain the 
(higher order in a s ) perturbative corrections, and 6*4(04), Cq(Oq) represent the leading non- 
perturbative corrections, of dimensions four and six, respectively [9]. They are dominated 
by the gluon condensate: C4(0 4 ) ~ -^-(a s G 2 ), and the four-quark condensate which, in the 
vacuum saturation approximation, is given by C 6 (0 6 ) ~ (l792/27N c )ir 3 a s (qq) 2 . Since the 
contribution of the condensates is negligible, and we have no new information to add, we 
simply accept the values quoted by BPR and JM/CPS in the remainder of this paper. 

To extract m u (s) +m ( i{s), one then needs to input the scale s = s , at which one assumes 
pQCD to have become valid and, second, experimental and/or model information on the 
hadronic spectral function (and hence its moments) below s . Having done so, one may then 
use either Eq. (3) or (4) to extract m u (s) + rrid(s), and from that, the MS combination of 
the masses at any desired scale p using the renormalization group running. Most sum rule 
analyses extract their estimates at ^fs > 1.7 GeV, and then run down to p — 1 GeV. We 
believe that it is unnecessary to introduce an extra uncertainty in the estimates by relying 
on pQCD over this interval where the running is large. For this reason our final comparisons 
are at p — 2 GeV. However, to preserve continuity with existing sum rule analyses masses 
quoted without any argument will always refer to the MS values at 1 GeV. 

The most up-to-date version of the above analysis was performed by Bijnens, Prades and 
deRafael[3] (BPR), whose treatment we will follow closely below. In this analysis, BPR have 
used the three-loop pQCD result of Refs.[10, 11] for the pseudoscalar two-point function, 
employing three active quark flavors with = 300 ± 150 MeV[12], and the values 

C 4 (0 4 ) = (0.08 ± 0.04) GeV 4 , (5) 
C 6 (0 6 ) = (0.04 ±0.03) GeV 6 . (6) 

for the non-perturbative, condensate contributions. For the hadronic spectral function on 
the interval (0, s) they include the pion pole, whose residue is known exactly in terms of 
f n and m n , and a 37r continuum contribution modulated by the n' and -n" resonances. The 
BPR ansatz is 

Phadronic(s) = ppole + F(s) P^pT^K 5 _ 9m l) ( 7 ) 

where the "37r continuum spectral function" P^p r j i (^) is obtained from the leading order, 
tree-level chiral perturbation theory (%PT) result for (0|<9 M A M |37r), and F is a modulating 
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Figure 1: The contour integral for the FESR's of the text. The "hadronic" integral from 
— > s, which includes contributions from the poles and cuts, is obtained using a model for the 
continuum portion of the spectral function, while the integral over the circle at sufficiently 
large s (s > sq) is done using the three-loop perturbative result. 



factor which accounts for the presence of the ir' and 7r" resonances. The form of F is taken 
to be a superposition of Breit-Wigner terms 

F( . A Wititts-Mt + iMiTiW* 

with £i = 1. 

There remain three unknowns at this point, the overall normalization parameter A, the 
relative strength and phase, £2, of the two resonances, and the value of sq . BPR find that, if 
they assume sq ~ 2-3 GeV 2 , duality can be satisfied for a number of values of [A, £2]- Their 
best solution uses the normalization A = 1 at threshold, and then fixes £2 by demanding 
duality between the hadronic ratio 

3 J°dttp 5 {t) 

n had .{s) = - js dtp5{ty (9) 
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and its pQCD counterpart 



Kqcd(s) 



l + R 1 (s) + 2^^' 



(10) 



over the interval between the two resonances, i.e. 2.2 < s < 3.2 GeV 2 . We have reproduced 
the results of BPR with their choice of resonance parameters (which differ slightly from those 
listed in their published version [13]) based on PDB94 [12], 



Their preferred solution (solution 2) is shown in Fig. 2. For s > so = 3 GeV 2 we also 
plot the perturbative spectral function (duality constraint) for m u + = 12 MeV, their 
extracted value of the quark mass. As is evident from the figure, the rise due to the 7r"(1800) 
is roughly consistent, both in magnitude and slope, with the perturbative ansatz. This is 
a consequence of tuning the normalization and relative phase of the second resonance, and 
leads to approximate duality over the range 2.2 GeV 2 < s < 3.5 GeV 2 . However, the fall-off 
of the spectral function on the far side of the 7r"(1800) resonance, in contrast to the rising 
pQCD solution, shows that, in order to preserve duality, further resonances and intermediate 
states are required to bolster the BPR ansatz beyond the 7r"(1800) peak. 

Note that, in the BPR analysis, the threshold behavior of the spectral function is not 
determined experimentally, but rather obtained from leading order ;\PT. To the extent 
that SU(2) x SU(2) ;\PT converges well at leading order, the choice A = 1 then ensures 
correct normalization of the spectral function near 3n threshold. However, in the spectral 
integral appearing in Eq. (3), which determines the light quark mass, the contribution of 
the near-threshold region is negligible compared to that from the vicinity of the resonance 
peaks. Correctly normalizing the spectral function in the resonance region is thus much more 
important than correctly normalizing it near threshold. We will show later, by considering 
an analogous example (the correlator of two vector currents), that the conventional threshold 
constraint, A — 1, almost certainly leads to a significant overestimate of the spectral function 
in the resonance region. 

To summarize, we will investigate the following aspects of the BPR solution: the un- 
certainty in the mass extraction produced by uncertainties in the 3-loop pQCD expression, 
the reliability of the overall normalization of the continuum contribution, and the sensitivity 
of the results to the value chosen for sq , the scale characterizing the onset of duality with 
pQCD. The same issues are also relevant to the extraction of m s using the Ward identity for 
the vector current. Our contention is that plausible systematic errors in each are such as to 
lower the estimates for light quark masses. 

3 Convergence of 2-point functions in pQCD 

The pseudoscalar 2-point function is known to three loops in pQCD [10, 11]. The main issue, 
in applying this expression to the problem at hand, is the question of convergence. If, for 
example, we write 1 + with i?j as defined in Eqs. (3) and (4), in the form (1 + xa s /ii + 
y(a s /Tr) 2 ), then the coefficients x,y show a geometrical growth, i.e. the growth for R x and 



M l = 1300 MeV, r x = 325 MeV; M 2 = 1770 MeV, T 2 = 310 MeV . 



(11) 
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i?2 is roughly the same and the average values are x ~ 6.5 and y m 46. As a result the 
0(a s ,a^) correction terms are 0.61 and 0.41 respectively at s = 3 GeV 2 where a s /n m 0.1. 
Since these are large, it is important to estimate the sum of the perturbation series. One 
plausible possibility is to represent the series by the Pade 1/(1 — 0.63), in which case the 
neglected terms would further increase the pQCD estimate by ~ 35% and consequently 
lower the BPR result for m u + by \A-35, i.e. from 12 — > 10.4 MeV. In fact, historically 
sum-rule estimates have decreased over time precisely because of the increase in the pQCD 
result. Part of the change has been due to the increase in the value of A.qq D an d part due to 
the large positive 3- loop contribution [3] . If this trend were to continue, and the unknown 
higher order terms were to continue to grow geometrically and contribute with the same sign 
(as is the case for the scalar channel discussed below) then, the extracted quark mass could 
be significantly lowered. 

The situation in the case of the scalar 2-point function analyzed by JM/CPS is somewhat 
better. A very careful analysis of the stability of the pQCD expressions and of the choice 
of the expansion parameter has been carried out by CPS [4] who include terms up to 
in the 2-point function and in the running of the coupling and mass. The pQCD result, 
after Borel transformation, has the expansion (1 + 4.8a s /7r + 22(a s /7r) 2 + 53(a s /7r) 3 ) [5]. 
Taking a s /ir m 0.1, as appropriate for u = 4 GeV 2 with Aq^ d = 380 MeV, we find that the 
difference between the pQCD series and a possible Pade representation, 1/(1 — 0.48), is only 
about 9%. This correction would lower the estimate of m s by ~ 5%, consistent with the 
estimate by JM [5]. 



4 Constraints on m u + rrid from the Posit ivity of p 5 (s) 

The fact that the spectral function ps(s) is positive definite above threshold allows us to 
place rigorous lower bounds on m u + rrid as a function of So [14]. A weak version of this 
bound (labeled "pole") is obtained by ignoring all parts of the spectral function except for 
the pion pole, whose contribution to the integral in Eq. (3) is 2/ 2 m^ (Eq. (4) produces a 
much less stringent bound and hence is not considered further). One then finds, assuming 
the validity of the input 3-loop pQCD result, 

where s is the upper limit of integration in Eq. (3). A stronger constraint (labeled "ratio") 
is obtained by noting that, for p 5 (t) > 0, 



IL dttp 5 (t) 



< s (13) 



i: th dt P5 (t) 

where sth denotes the 37r threshold value. The bound is saturated when the entire spectral 
strength is concentrated as a delta function at s. If s in Eq. (13) is assumed to be in the 
dual region, this turns out to place considerably stronger constraints on m u (s) + m d (s). To 
see this, note that the LHS of the inequality in Eq. (13) is, using Eqs. (3) and (4), 

^K(s) + m d {s)f f {l + R 2 (s) - l^el} - 2/M 



87I- 2 



K(s) + m d {s)Y f {l + Ri(s) + 2^i> } - 2/>4 • 



(14) 
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From the expression (14) we see that, if we start with a large value of m u (s) + m^(s) and 
begin to lower it, keeping s fixed, the inequality (13) will be violated before we reach the 
value of m u (s) + m^(s) corresponding to the pion pole saturation of the spectral function 
at which point the denominator in (14) vanishes). Thus the inequality (13) provides a more 
stringent (larger) lower bound on the extracted quark mass. This is illustrated in Fig. 3 
where the dependence of (m u + majnUn on s for both of the above constraints is shown. 

The "ratio" curve shows that if one assumes so — 2.5 GeV 2 , as in the BPR analysis, then 
n^u + m d > 10 MeV. The fact that the BPR result for the mass extraction, m u + m d ~ 12 
MeV, is close to this lower bound is a reflection of the fact that the spectral strength is 
concentrated in the region close to the assumed onset of duality. Such a feature is, in 
fact, rather natural since Sq is chosen to coincide with the 7r"(1800) peak. However, if sq is 
considerably larger than 3 GeV 2 (to alleviate the problem of large 0(a s , a 2 ) corrections to Ri 
at s ~ 3 GeV 2 discussed above) then considerably smaller masses are allowed by the "ratio" 
constraint, as is evident from the figure. Furthermore, one would, in fact, expect masses not 
much greater than the "ratio" bound to be favored in all cases where spectral functions are 
characterized by resonance modulation of a rising continuum phase space background and 
have their spectral strength concentrated in the region near s - 

These bounds make it clear that the value of the quark mass extracted from FESR 
will tend to be very strongly correlated with assumptions about the appropriate value of 
s . In addition, it will, of course, depend on the details of the hadronic spectral function 
from 37r threshold up to So, which are, at present, not experimentally determined. Since the 
perturbative p§{t) is known up to the overall normalization, which is given by the quark mass, 
one test of the validity of the phenomenological ansatz for the hadronic spectral function 
would be to show that the results for the quark mass remained stable under variations of the 
upper limit of integration, s, in Eqs. (3) and (4). This test, however, is meaningful only if one 
already knows that the values of s being employed are greater than sq. Unfortunately, the 
lack of experimental information on the hadronic ps(s) precludes the possibility of making 
such a test. In the next Section we construct plausible spectral functions, all satisfying 
duality, corresponding to a range of possible values for sq lying between 3 and 10 GeV 2 , by 
including higher resonances in the 37r channel. These models illustrate how, in the absence 
of experimental information, the uncertainty in m u + might be as large as a factor of 2 
if considerably higher values of s are chosen. 

5 Plausible Spectral Functions ps(s) in the Dual Region 

The assumption of duality places constraints on the form of the spectral function, ps(s). 
Below s = so, these constraints amount only to the determination of certain moments of 
the spectral function on the interval (s t h, s ), and hence are not particularly strong. In fact, 
as we illustrate below, the constraints of Eqs. (3) and (4) allow considerable freedom in the 
choice of ps(s) for s < s . For s > Sq, in contrast, duality determines the "average" Ps(s), 
i.e. averaged over some suitable region of s. This average value is given by the perturbative 
P5(s) , which can be obtained straightforwardly by differentiating the RHS of either Eq. (3) or 
(4). We evaluate these derivatives numerically using either of the two forms, which of course 
give consistent results. Even if one eliminates the running masses by matching the ratios 
of Eqs. (9) and (10) for s > s , it is easy to show that the resulting equation completely 
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determines the perturbative Ps(s), up to an overall multiplicative factor, for all s > s Q . 
The result of the duality constraints, in either form, is that ps(s) must be a monotonically 
increasing function of s, for s in the duality region. Numerically we find that this function 
is approximately linear as illustrated in Fig. 2. 

The hadronic spectral function in this channel is not known experimentally. It receives 
contributions not only from the pion and its resonances, but also from the resonant and non- 
resonant portions of the 3tt, 5tt, 7tt, KK71, . . ., NN, . . . intermediate states. Experimentally, 
only the 7r'(1300) and 7r"(1800) have been observed as distinct resonances [15]. Even so, 
their decay constants are not known experimentally, and hence the normalization of their 
contributions to the spectral function have to be treated as free parameters. The number of 
multi-particle intermediate states one has to consider, moreover, grows with s , as does the 
problem of separating their resonant and non-resonant portions. From dimensional argu- 
ments, the contribution of these various intermediate states will grow linearly at sufficiently 
large s. In the region of resonances, the resonances will modulate the cut contribution, 
and the hadronic spectral function is expected to match the pQCD behavior only after an 
average over some interval of s. This averaging is crucial if the resonances are narrow and 
isolated. Alternately, if the widths of subsequent resonances, for example n"' and it"", be- 
come much greater than the resonance separation, then the overlap of resonances can provide 
the monotonically rising behavior required by duality, and averaging is not crucial. 

We illustrate these points by constructing ansatze for the hadronic spectral function 
which are of the form used by BPR, i.e., involving resonance modulation of the continuum 
37r background. To explore values of s as large as 10 GeV 2 with the ansatz above, we 
include pseudoscalar resonances with masses as large as ~ vTU GeV. For the first two such 
resonances, the 7r'(1300) and 7r"(1800), we use the 1996 Particle Data Group values for the 
masses and widths. For the remaining two resonances, the n'" and n"", expected in this 
range, we are guided by model predictions. The n'" resonance is typically expected to lie 
around 2400 MeV in models constrained by the lower part of the meson spectrum[16]. In 
addition, the 3 Po model[17], which has proven reasonably successful in estimating decay 
widths[18], predicts a width for the 7r'"(2400) between 700 - 1900 MeV[16, 19] depending 
on how the relativistic effects are treated. The approach leading to 700 MeV gives 300 MeV 
for the width of the 7r"(1800) which is larger than the experimental value of 212(37) MeV. 
We, therefore, assume the lower limit 700 MeV for the width in this study, even though this 
may be an overestimate. Similarly, we assume that it"" lies at 3150 MeV with a width of 
900 MeV. In short, we choose 

Mi = 1300 MeV, F x = 325 MeV; M 2 = 1800 MeV, T 2 = 212 MeV; 

M 3 = 2400 MeV, T 3 = 700 MeV; M 4 = 3150 MeV, T 4 = 900 MeV. (15) 

The decay constants of all of these resonances are unknown, and will therefore be treated as 
free parameters. The limitations of such a truncated spectral function are obvious; however, 
it should be noted that, because we have allowed ourselves some phenomenological freedom 
in treating the strengths and widths of the last two resonances, our ansatze for the spectral 
function can also be thought of as providing an approximate means of representing a combi- 
nation of resonant and non-resonant effects. Our aim is, in any case, to simply demonstrate 
how the piling up of resonances can give the pQCD behavior, and the nature of plausible 
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Table 1: The parameters used to generate the spectral functions shown in Fig. 4. The 
normalization at threshold A and the relative weights q assigned to the resonances are 
defined in Eq. (17). 







m u + m d ( MeV) 


Cl 


c 2 


c 3 


c 4 




Case 1 


3.0 


12.0 


1 


-0.23 + 0.65i 






1.0 


Case 2 


5.7 


9.0 


1 


1.0 


2.3 




1.0 


Case 3 


8.0 


8.0 


1 


1.2 


5.0 


6.5 


0.7 


Case 4 


10.0 


6.0 


1 


0.8 


2.0 


3.68 


0.5 



spectral functions for which the "extracted" quark mass is, as for the BPR case, rather close 
to the value given by the "ratio" bound. 

For the resonance modulated spectral function we adopt, following BPR, the ansatz 



Ps(s) 



F(s) p^ FT {s) 

where p 37 pp(s) is the spectral function corresponding to the leading order, tree-level 
result for (0|c^A^|37r), and 

J2i CiMiTi/[(s - M 2 ) 2 + M?rf] 



(16) 



F(s) 



A 



EicMFi/Ks^-Miy + Mir 2 } 



(17) 



The sum in Eq. (17) runs over the appropriate number of resonances, depending on s as 
described below, with relative strengths q. The parameter A is the overall normalization of 
the resonance contribution to the continuum part of the spectral function at 37r threshold. 
We have taken the q to be real, in order to simplify the task of searching for suitable spectral 
functions, whereas BPR, who use a slightly different form for F, as given in Eq. (8), with 
just the first two resonances, allow the relative strength of the two resonances to be complex. 

We display a series of spectral functions in Fig. 4, all satisfying duality and constructed 
by employing up to four resonances in the ansatz above. The values for A, {q}, Sq and 
m u + m d used in the construction are given in Table 1. As one can see from the figure, 
there exist perfectly plausible spectral functions corresponding to m u + m d = 12, 9, 8, and 
6 MeV. The first case (s = 3 GeV 2 , m u + m d =12 MeV) is the BPR solution discussed 
before. The second case (s = 5.7 GeV 2 , m u + m d = 9 MeV) corresponds to including three 
resonances and matching to the duality solution at the top of the third resonance. The 
assumption here is that the third and higher resonances merge to produce the dual solution 
above this point. The matching in the third case (s = 8 GeV 2 , m u + m d = 8 MeV) is at 
the beginning of the rise of the fourth resonance, while in the fourth case (s = 10 GeV 2 , 
m u + md = 6 MeV) we match at the top of the fourth resonance. In cases where we match 
at the peak of a resonance, the dual region actually appears to begin somewhat below the 
input value of Sq . This is because the slope of the rising side of the last resonance tends to 
match reasonably well the slope of the pQCD version of P5(s). 

These spectral functions are, by construction, perfectly dual for s > s . Duality also 
requires the low-energy (s < s ) part of ps(s) to have the correct moments to satisfy Eqs. (3), 
(4). However, the constraint of duality does not lead to a unique solution. Experimental 
data (decay constants) are needed to fix the overall normalization A and the relative weights 
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Table 2: The parameters used generate the plots shown in Fig. 5. The values of s and 
m u + m d have been fixed to s = 10 GeV 2 and m u + m d = 6 MeV respectively in each of 
the three cases. 





Cl 


c 2 


c 3 


c 4 


A 


Case 1 


1 


0.8 


2.0 


3.680 


0.5 


Case 2 


1 


0.6 


2.0 


5.160 


0.4 


Case 3 


1 


1.0 


4.0 


11.75 


0.3 



Cj. We illustrate this point in Fig. 5 by constructing three spectral functions that differ for 
s < s . In all three cases s and the input m u + m d in the pQCD expression are fixed 
to be the same, while the values of parameters A and q are as defined in Table 2. The 
corresponding output values for m u + m d are shown in Fig. 6. As expected they converge 
to the input value in the dual region. 

There are three features of these ansatze that should be noted. First, the value of 
Tfin + m d decreases with s in a manner very similar to the "ratio" bound. This is because 
in each case the spectral function is stacked up towards so- Second, we find that, to produce 
spectral functions corresponding to values of m u + m d only a few MeV above the "ratio" 
bound, the threshold normalization parameter, A, has to be decreased with increasing sq. 
In the next section we will show that values of A significantly smaller than 1 are, in fact, 
to be expected, based on a consideration of the analogous vector current correlator, for 
which the normalization in the resonance region is known experimentally. Third, the q 
are large. It is not clear, a priori, if this should be considered unreasonable or not. For 
example, in the narrow width approximation the q would scale as ~ {ffMf) / (f%M*) and 
thus have an explicit dependence on Mf. (The BPR model spectral functions, being even 
larger than ours, of course, correspond to even larger 7r'(1300) and 7r"(1800) decay constants.) 
Moreover, by leaving the normalizations as free parameters we are potentially incorporating 
other non-resonant background effects. Ultimately, this issue can only be resolved by appeal 
to experimental data which, unfortunately, is not available at present. 

The bottom line of the above discussion is that since both the correct value for the location 
of the onset of duality with pQCD, and the correct form of the hadronic spectral function, 
are at present unknown, the value of m u + m d extracted using FESR can easily vary by a 
factor of 2. As we have pointed out, using s ~ 3 GeV 2 leads to a perturbation series in which 
the a s and ct 2 s terms are large. As soon as one allows significantly larger values of Sq, in order 
to alleviate this problem, however, considerably smaller values of the extracted quark mass 
are possible. We will now, furthermore, argue that the conventional method of normalizing 
the continuum part of the spectral function tends to produce significant overestimates of 
the resonance contributions, and hence also significant overestimates of the extracted quark 
masses. 
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6 Normalization of the resonance-modulated spectral 
function 



The assumption that a given spectral function may be written as a Breit-Wigner resonance 
modulation of a continuum phase space factor, as exemplified in Eqs. (7) and (8), is valid 
in the vicinity of a narrow resonance. We will assume that this ansatz, used by both BPR 
and JM/CPS (see section 7) in the pseudoscalar and scalar channels respectively, is a good 
approximation. This fixes the general form of the spectral function but, of course, does 
not fix the magnitude in the resonance region, since the relevant pseudoscalar and scalar 
decay constants (with the exception of /„-) are not known experimentally. Both BPR and 
JM/CPS deal with this problem by assuming that resonance dominance of the relevant spec- 
tral function continues to hold all the way down to continuum threshold. Thus, for example, 
the overall scale of the BPR and JM/CPS ansatze for the continuum part of the spectral 
function is obtained by choosing A — 1, i.e., by assuming that the tails of the resonances 
reproduce the full threshold spectral function. BPR, in the absence of experimental data, 
use the tree-level ;\PT expression for the spectral function in the threshold region. The 
JM/CPS treatment differs only in that they normalize the sum-of-resonances ansatz at Kir 
threshold using experimental data (the scalar form factor at threshold is computed using the 
Omnes representation with experimental Kn phase shifts as input.) 

The second key point in the JM/CPS ansatz for spectral function is the assumption that 
one can take the "standard" s-wave s-dependent widths for the resonance contributions. This 
assumes that the effective coupling of the strange scalar resonances to Kn is momentum- 
independent over the whole kinematic range relevant to the spectral integral. We will now 
show that the combination of this assumption and of the resonance-saturation-at-threshold 
can fail badly by studying its exact analogue in the isovector vector channel. In fact, in the 
vector channel, the analogous set of assumptions produces a significant overestimate of the 
spectral strength in the region of the resonance peak. 

Consider, therefore, the vector correlator 



n 33 (g 2 ) = (qi^Qu - q 2 g»v) n 33 (g 2 ) 

= i J d A xe^(0\T{Vr(x)V»(0)}\0) , 



(18) 



where is the 1 = 1 vector current. In the narrow width approximation, the p contributions 
to the spectral function of n 33 , p 33 , at the p peak, is known in terms of the p decay constant, 
F p = 154 MeV, 



F 2 

p 



nT p m p 



0.0654 . 



(19) 



Let us now apply the analogue of the BPR and JM/CPS ansatze to the vector channel, by 
assuming the (trial) spectral function to be given by 



48tt 2 



Ami 



L C(K) 



(20) 
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where the quantity in the first set of square brackets is the leading order x?T expression for 
the spectral function [20], and 



sm p T p / 


m 2 , — 4m? 


3/2 


{m? p - s 


) 2 + I»2m2 



(21) 



which follows from employing the p-wave s-dependent width 



m p Y p 



(m 2 p - 4ml) 



3/2 



1 - 



4ml 



3/2 



s . 



(22) 



We have chosen this form of the width in analogy to the "standard" s-wave s-dependent 
width of JM/CPS. This ansatz assumes that the effective coupling of the p to -kit has the 
minimal form g p7T7T p p (7r + <9 M 7r~ — 7i^dii + ) with g pn7T independent of momentum over the rele- 
vant kinematic range. The threshold factor [1 — 4m 2 r /s] 3 ^ 2 in the numerator of Eq. (22) has 
been separated out explicitly in writing Eq. (20). The ansatz Eq. (20) then implies 



/47'K) 



m 



2 _ 4 m 2 



7/2 



1927T 2 



m 3 p T 2 p ml 



0.27 , 



(23) 



a factor of 4.1 too large. Had we instead used the normalization given by the full next-to- 
leading order xPT expression [20], (which matches well to experimental data near threshold), 



48tt 2 



3/2 



0(s - 4ml) 



' | 4LS(/x) 8 | 

f2 ~ l ~ 

J 7T 



(24) 



the peak height would be further increased by a factor of 1.28, the correction being dom- 
inated, for /i ~ m p , by the term in the square brackets in Eq. (24) involving the 0(q 4 ) 
renormalized low-energy constant (LEC), L r 9 , 



4U 9 {m p ){4ml) 



= 0.24 , 



(25) 



where we have used L r 9 {m p 



0.0069(2) [21], and + • • • refers to loop contributions whose 

form is not important in what follows. Note that, since it is the next-to-leading order 

expression, Eq. (24), which matches experimental data, it is this latter normalization which 

corresponds to the JM/CPS treatment of the scalar channel. The analogue of the JM/CPS 

ansatz, in the case of the vector correlator, thus overestimates the spectral function at the p 

peak by a factor of 5.1. 

The source of this problem is not difficult to identify and, in fact, turns out to be that 

the crucial assumption that the spectral function can be taken to be completely resonance 

dominated, even near threshold, is incorrect. This is most easily seen from the perspective 

of xPT. Indeed, it is known that, when one eliminates resonance degrees of freedom from 

a general, extended effective Lagrangian, producing in the process the usual effective chiral 
yPT 

Lagrangian, , relevant to the low-lying Goldstone boson degrees of freedom alone, the 
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effect of the resonances present in the original theory is to produce contributions to the 

yPT 

LEC's appearing in C£ff [22, 23]. There are two important observations about the nature 

of these contributions which are of relevance to the present discussion. The first is that 

yPT 

the resonances do not contribute to the lowest-order (0(q 2 )) LEC's of C^ff ', instead, the 
leading (in the chiral expansion) contributions are to the 0(q 4 ) LEC's, L r k (p) (where p is 
the yPT renormalization scale, and we adhere throughout to the notation of Gasser and 
Leutwyler[24]). The second is the phenomenological observation that, if one takes p ~ m p , 
the resonance contributions essentially saturate the L k (p)[22, 23] (see, for example, Table 
2.1 of Ref. [21], for a comparison with recent experimental determinations of the LEC's). 
An immediate consequence of the first observation is that the correct normalization for the 
resonance contributions to quantities like P33, or (0|<9^A M |37r), near threshold, cannot be that 
coming from the tree-level (0(q 2 )) yPT contributions, since such contributions are associated 
with the Goldstone boson degrees of freedom alone, and contain no resonance contributions 
whatsoever. Similarly, normalizing to the full threshold value, as obtained, for example, 
from experiment, would also be incorrect, since this full value necessarily contains both tree- 
level and leading non-analytic contributions, neither of which can be associated with the 
resonance degrees of freedom, in addition to the 0(q A ) LEC contributions which do contain 
resonance contributions. Fortunately, the second observation provides us with an obvious 
alternative for normalizing resonance contributions near threshold. We propose, therefore, 
to accept the phenomenological observation above as a general one and identify resonance 
effects in near-threshold observables with those contributions to the 1-loop expressions for 
these observables involving the appropriate 0(q 4 ) LEC's, {L r k }, evaluated at a scale p ~ m p . 
Such an identification, however, requires that the LEC is dominated by the appropriate 
resonance, as is the case for the vector (Lg) and scalar (L 5 ) channels, but not for the pseudo- 
scalar channel. This prescription, like that of BPR and JM/CPS, represents a means of using 
information solely from the near-threshold region (in this case, obtainable from a knowledge 
of the chiral expansion of the spectral function) to normalize the spectral function in the 
resonance region. However, we will show below that, in contrast to the analogue of the 
BPR and JM/CPS ansatze which was in error by a factor of ~ 5 at the p peak, the new 
prescription normalizes the peak accurate to within a few percent. Based on the success 
of the prescription in this channel, we will then apply it to a re-analysis of the JM/CPS 
extraction of m s involving the correlator of the divergences of the vector current. 

Let us return, then, to the spectral function p 33 . According to the discussion above, the 
p meson contributions to P33, near threshold, can be obtained by taking just that term in 
Eq. (24) proportional to Lg, evaluated at a scale p ~ m p . The only change in the above 
analysis is then a re-scaling of p^™ 1 in Eq. (23) by a factor of 0.24, the value of the 0(q A ) 
LEC contribution in Eq. (24) at p = m p . This leads to a prediction for the spectral function 
at the p peak of 

P L 33 EC H) = 0.067 (26) 

in good agreement with the experimental value given in Eq. (19). 

Let us stress that the precise numerical aspects of the prescription above, namely the sup- 
position that the normalization at resonance peak of resonance contributions to the hadronic 
spectral function can be obtained by evaluating the relevant 0(q A ) LEC contributions ap- 
pearing in near-threshold yPT expressions, at a scale p ~ m p , is one that is purely phe- 
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nomenologically motivated [22, 23]. While highly successful in the case of the vector channel, 
it has not been tested outside this channel. The fact that resonance contributions begin only 
at 0(q 4 ) in the chiral expansion, and hence that resonances do not contribute to either 
lowest-order tree-level or leading non-analytic terms in the x?T expansions of the relevant 
spectral functions, however, clearly indicates, independent of the numerical reliability of this 
prescription, the unsuitability of normalizing the resonance peaks by associating the full 
%PT or experimental values near threshold with resonance effects. Moreover, so long as 
the spectral functions of interest have even reasonably normal chiral expansions, with the 
dominant contributions near threshold coming from the lowest order tree-level contributions, 
we can conclude that the standard method of normalization will produce values for these 
spectral functions in the resonance region that are overestimated by a significant numerical 
factor. 

At this stage we should also mention that Stern and collaborators have suggested that 
the normalization at threshold could actually be much larger than that given by leading 
order xPT, as is expected in "generalized xPT"[25]. They then argue that, in that case, 
the quark masses would be even larger. Our observations are also relevant in this case: we 
again stress that, since the sum rules we consider are dominated by the resonance region, 
threshold normalization will only provide useful input if one can disentangle the contributions 
to threshold amplitudes associated with resonances from those associated with the Goldstone 
boson degrees of freedom. 

7 Re-analysis of the JM/CPS extraction of m s 

In this section we will employ the prescription proposed above to a re- analysis of the JM/CPS 
extractions of m s [4, 5]. Such a re-analysis is possible in this case because the 1-loop %PT 
expression for the relevant scalar form factor is known [26]. To introduce notation, we 
briefly review the analysis of Ref. [5] (that of Ref. [4, 27] is similar and need not be discussed 
separately). These analyses involve a standard QCD sum rule treatment of the correlation 
function 



where V^(x) is the strangeness-changing vector current and S(x) the corresponding strangeness- 
changing scalar current. The correlator of scalar currents is evaluated using the opera- 
tor product expansion (OPE). All terms on this side of the sum rule are proportional to 
{m s — m u ) 2 , and the full pQCD result is known for the predominant contribution \l/ [4]. 
The hadronic spectral function in the phenomenological side is again taken to be a sum-of- 
resonances modulation of the spectral function relevant to the Kn intermediate state near 
threshold. 

JM/CPS write the Kir contribution to the physical spectral function as 




{m s — m. 



u 



) 2 i I d A x e iq - x (0\T{S(x)S\0)}\0), 



(27) 



= ^n 9 ^ - s +)\A s ~ s +)( s - s -) \ d ( s )\ 2 



(28) 
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where s± = {tuk ± inn) 2 , and d(s) is the strangeness-changing scalar form factor, measured 
in Ki 3 for mj < s < s_, 

d(a) = (m 2 , - mj)/ (a) = (rr& - m 2 )/+( S ) + s f_(s) (29) 
with /±(s) the usual form factors defined by 

^{p ! )\- s1 mk + {p)) = [(jZ+pUW + (p - pV-(*)] • (30) 

In their analysis, JM/CPS employ the following resonance-modulation ansatz for the spectral 
function: 



Phadronic(s) = T^^V '( S ~ S +)( S ~ S -) K S +) P F ( S ) ( 31 ) 

where 



y c BW (s) 

F (s) = IT") \ (32) 

with 



r .BW( \ _ fn m n^n /qo\ 

In Eqs. (31)- (33), s + is the continuum Kn threshold, and f n , m n , and T n are the decay 
constant, mass and width of the nth scalar resonance, T n (s) being the usual s-dependent 
width given in [5]. The s-dependence of the width factor occurring in the numerator of the 
Breit-Wigner resonance forms has already been factored out explicitly in writing Eq. (31). 
The sum in Eq. (32) is taken to run over two resonances (the _K"q(1430) and Xq(1950)), and 
the duality point, s , of QCD sum rules (describing the point beyond which the physical 
spectral function is to be modeled by its perturbative expression) is fixed by a stability 
analysis. Note that the normalization procedure above assumes that the physical spectral 
function is completely saturated by resonance contributions near threshold. The threshold 
value of the scalar form factor, d(s+) = 0.33 ± 0.02 GeV 2 , is obtained using Omnes repre- 
sentation with experimental Kn phase shifts as input. This result is, moreover, shown to be 
consistent with that of %PT to l-loop, which can be obtained from the expression for fo(s) 
given by Gasser and Leutwyler[26] (g?q PT (s + ) = 0.35 GeV 2 ). Lastly, the master equation 
used for extracting m s is [5] 

u ^ope = / e s/u p h adronic ds + e s/u ppQCD ds (34) 

JO J so 

where both ^>" OPE and p p qcd are proportional to (m s — m u ) 2 . 

The first of the three issues raised by us, namely the reliability of pQCD has already 
been discussed in Section 3. We agree with JM/CPS that in this channel the effect of the 
neglected a A s and higher contributions could, at best, lower estimates of m s by ~ 5%. The 
remaining two issues, the value of s and the normalization of the hadronic spectral function 
are far more serious, as we now explain. 

To elucidate the role of s in the JM/CPS analysis we plot, in Fig. 7, both the model 
JM hadronic spectral function (for s < sq) and the pQCD version of the spectral function 
(for s > sq). We have used the JM values corresponding to the preferred solution, i.e., 
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s = 6.0 GeV 2 , A ( q ] cd = 380 MeV and m s = 189 MeV. The plot shows very clearly that the 
ansatz for phadronic is, at best, valid only for s < 4.0 GeV 2 . Furthermore, as evident from 
Eqs. (31)-(33), phadronic goes to a constant at large s, whereas p p qcd grows linearly (with 
logarithmic corrections). For this reason there is a large discontinuity between phadronic and 
Ppqcd even for s as low as 4 GeV 2 . The only way that Phadronic constructed from the Kit 
channel can satisfy duality is if there is a piling up of higher resonances, and these have 
to have large amplitudes (as we illustrated in section 5 for the pseudoscalar channel). We 
contend that so should only be chosen in the range where phadronic is known reliably. However, 
for s < 4.0 GeV 2 , and using the JM/CPS ansatz for phadronic, we have not been able to find 
a result for m s that is stable under variations of the Borel parameter u. It was precisely 
this lack of stability that forced JM/CPS to choose a larger s . Such a choice, we contend, 
is not reasonable as Phadronic <S P p qcd over the range 3 < s < 6 GeV 2 , i.e., duality is badly 
violated over this whole range. 

Lastly, we turn to the quantity, d(s + ), which sets the overall normalization of the reso- 
nance contributions in Eq. (31). This quantity is crucial in the JM/CPS analysis since, as 
noted by JM, the extracted value of m s scales directly with d(s+). The problem is that, just 
as for the light quark case, the spectral integral appearing on the phenomenological side of 
the sum rule is dominated, not by the near-threshold region, but by resonance contributions. 
The ansatz (31)-(33) for the spectral function, however, is designed only to produce the 
correct overall normalization at Kir threshold. From our discussion above of the analogous 
treatment of the vector current correlator, it is clear that such an ansatz will overestimate 
the resonance contributions near threshold, and hence almost certainly significantly overes- 
timate the spectral function in the resonance region. To correct this problem we need to 
properly rescale the JM/CPS ansatz at threshold. We do so on the basis of the proposal 
above, i.e., we assume that in the scalar channel, just as in the vector channel, the 0(q A ) 
LEC's, evaluated at a scale p ~ m p , give the correct normalization of the scalar resonance 
contributions at threshold. It is easy to implement this revised normalization of phadronic 
because, not only is the 1-loop expression for d(s) known[26], but, in addition, Jamin 
and Miinz have demonstrated explicitly the accuracy of this expression for d(s + )[5]. 

Let us write the 1-loop expression for d(s + ) in the form 

tree(s+) ~\~ < ^res(^+) P) + di oop {s + ,p) (35) 

where d tre e(s+) is the leading, 0(q 2 ) tree-level contribution, d res (s + ,p) contains the 0(q A ) 
LEC contributions, and di oop (s + ,p) contains the contributions associated with 1-loop graphs 
generated from the 0(q 2 ) part of the effective chiral Lagrangian. The latter two terms 
are separately scale-dependent. According to the prescription introduced above, resonance 
contributions to d(s + ) are to be identified with d res (s + , m p ). Resonance contributions to 
|(i(s + )| 2 , consistent to 1-loop order, are thus given by 

\d(s+) \ 2 res ~ 2 d tree (s + ) d res (s + , m p ) . (36) 

Using [26] 

d tree {s+) = {m\ - ml) = 0.22 GeV 2 
d res (s + , m p ) = As + (m 2 K - ml)L r 5 (m p )/f 2 , (37) 
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with L r 5 (m p ) = 0.0014±0.0005, we find that \d(s + )\ 2 res ~ 0.23|rf(s + )| 2 . With no other changes 
to the JM/ CPS analysis than the corresponding re-scaling of the continuum spectral function, 
the value of m s would thus be lowered by almost exactly a factor of 2. However, as discussed 
above, there are problems of consistency with using the JM/CPS ansatz for the spectral 
function with values of s as large as 6 GeV 2 . 

In light of the above corrections, the question before us is whether it is possible to get a 
stable estimate of m s by repeating the JM/CPS analysis with sq m 4 GeV 2 and an overall 

(3) 

normalization of Phadronic oi A ^ 0.25. To answer this question we have varied A.q CD in the 
range 200 — 450 MeV, the relative strength /2//1 of the two Breit-Wigner resonances in the 
modulating factor over 0.2 — 1, and A over the range 0.2 — I. Despite this, we have failed 
to find a solution that is stable under variations in the Borel scale u. The cause of this 
failure is the ansatz for phadronic, and the small range of s over which it can be evaluated. 
It is our contention that reliable results for m s using sum-rules can only be obtained if 
Phadronic is determined to high precision over a sufficiently large range of scales, say from s t h 
to 8 GeV 2 . If s is "small" then limitations of operator product expansion, convergence of 
perturbation theory at small s, and details of resonances contributions make it difficult to 
test the reliability of the results. 

For completeness we should also mention the alternate JM determination of m s via an 
analysis of the analogous strangeness-changing axial correlator. Their result in this case is 
m s (MS, 2 GeV) = 91 MeV, significantly smaller than that obtained from the scalar channel 
via the treatment of the vector current correlator. They, however, consider this analysis 
incomplete because it employs, for the normalization of the continuum spectral function at 
threshold, the leading-order, tree-level xPT result. They contend, based on the expectation 
that the full normalization will, as in the scalar channel, significantly exceed that given by 
tree-level xPT (d(s + ) = 1.5d tree (s + ) for the scalar channel), that the true normalization will 
likely be significantly larger . If true, this would mean that m s would be correspondingly 
increased. They thus expect their two analyses to become consistent once they employ a 
normalization at threshold corresponding to the 1-loop expression for the continuum spec- 
tral function in the pseudoscalar channel. Our contention is that, in fact, the "correct" 
normalization is given, not by the full threshold spectral function, but rather by the appro- 
priate 0(q 4 ) LEC contributions to the 1-loop result, and that it should hence be significantly 
sma//er than that corresponding to the tree-level result. Further progress on this issue, and 
that of the consistency of the two different extractions for m s , will be possible only once 
the 1-loop expression for < Old^A^lK-mr > is known [28]. A re- analysis of the BPR FESR 
treatment of m u + md is similarly stymied by the absence of 1-loop expressions for the matrix 
elements < 0|9 M AW|37r >, and by the lack of association of the Lj involved with just the 
pseudoscalar resonances. 

In the past, of course, the agreement of the ratio r ~ 2(180)/12 = 30 obtained from the 
different sum rule analyses with that (24.4 ± 1.5) obtained from xPT[29] has been taken to 
provide a posteriori support for the validity of the sum rule treatments. Our contention is 
that a self-consistent sum-rule analysis would yield estimates of both m s and m u + m d that 
are lower by a factor of ~ 2, thus maintaining the consistency with the value of r. 
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8 Conclusions 



We have shown that the ability to make reliable extractions of m s and m u + m d using sum 
rule analyses rests on three key features of these analyses: the degree of reliability of pQCD, 
a knowledge of the scale, s , at which quark-hadron duality becomes valid, and an ability to 
construct hadronic spectral functions which are correctly normalized in the resonance region, 
even in channels where experimental data on the relevant decay constants is not available. 

We find that, in the relevant pQCD expressions, the a s , a 2 , . . . corrections are large both 
in the scalar and pseudoscalar channels. Including reasonable estimates for the unknown 
higher order terms lowers the sum rule estimates of quark masses. The largest effect is in 
the extraction of m u + md, which we estimate would be lowered by ~ 20% compared to the 
value quoted by BPR [3] . The correction in the case of m s extracted from the scalar channel 
is roughly 5%, and this has been accounted for by JM/CPS. 

Second, the estimates obtained for the quark masses are potentially very sensitive to the 
choice of s . We have illustrated this through an analysis of rigorous lower bounds and the 
use of a variety of plausible spectral functions in the case of m u + m^. Current sum rule 
analyses are forced to choose low values of s due to lack of experimental information. The 
FESR extraction of m u + m<z, for example, is based on rather low values of so < 3 GeV 2 , so 
no tests of the stability of the estimates under variations of s can be made. In the case of 
JM/CPS analysis of m S) the value chosen, s = 6.0 GeV 2 , is artificially large. This choice 
arises from an attempt to achieve stability of the Borel transformed sum rule with respect 
to the Borel parameter u. Since, however, the phenomenological ansatz for the spectral 
function breaks down for s > 4.0 GeV 2 , it is clear that such a choice of s is not physical. 
For reasonable choices of s we are also not able to find a solution that is stable with respect 
to variations in u. We, therefore, conclude that no reliable estimates of m s can be made 
unless Phadronic is known accurately over a significantly larger range of s. 

Third, we have shown that the method employed in previous analyses for fixing the over- 
all normalization of the resonance-modulated model spectral functions leads to significant 
overestimates of the continuum contributions to the relevant spectral integrals and hence to 
significant overestimates of the quark masses. The source of this problem is the fact that 
normalizing the resonance-modulated ansatz (see Eqs. (31)- (33)) to either the experimental 
value or to the %PT value for the spectral function in the near-threshold region results in the 
inclusion of near-threshold contributions of the Goldstone-boson degrees of freedom in addi- 
tion to the desired resonance contributions. Overestimating the size of the resonance tail in 
this manner, of course, leads to a corresponding overestimate of the resonance contributions 
at resonance peak. Since it is the resonance peak region, and not the threshold region, which 
dominates the phenomenological side of the sum rules, the conventional procedure produces 
significant overestimates of the quark masses. In the case of the vector current correlator, 
where the normalization of the spectral function at the p peak is known experimentally, we 
have shown that the magnitude of this overestimate is large: the conventional method of 
normalization produces a spectral function which, at the p peak, is a factor 4.1 — 5.1 larger 
than that given by experiment. We have explained, based on an understanding of the man- 
ner in which resonance effects manifest themselves in x?T why the conventional method of 
normalization cannot be correct, and have proposed an alternate phenomenological prescrip- 
tion for normalizing the spectral function, designed to provide estimates which are reliable 
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not so much in the threshold region as in the resonance region relevant to the sum rule quark 
mass extractions. We verify that this prescription reproduces the experimental result for the 
vector (p) channel. This method is straightforward to apply to the scalar channel as the 
1-loop (0(g 4 )) %PT corrections are known, and the revised estimate for the normalization 
could reduce the estimate of m s by as much as a factor of ~ 2 over the values found in 
previous analyses. We argue that a similar overestimate of the normalization will exist in 
the pseudoscalar channels, though we are unable to estimate its magnitude at present. 

The bottom line is that unless the hadronic spectral function is known accurately over a 
large range of scales, say up to s = 8 GeV 2 , reliable extraction of quark masses from sum- 
rules considered is not possible. Even though the lattice QCD estimates have their share of 
statistical and systematic errors [2], we claim that at present they represent the most reliable 
means of estimating the quark masses. Our estimates of the systematic errors in sum rule 
analysis suggest that revised sum rules estimates could easily be smaller by a factor of two, 
in which case these would be consistent with the small values obtained from lattice QCD. 
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Figure 2: The hadronic spectral function assuming that the dual region begins at the ir" 
resonance. The location, widths, and normalization are the same as in the solution found by 
BPR. The dashed line is the continuation of the BPR solution for s > sq. For s > sq we also 
show the spectral function required to satisfy perturbative duality for m u + m d = 12 MeV. 
The two ansatze are joined smoothly by choosing s = 3.0 GeV 2 . s is in GeV 2 , and p 5 in 
GeV 4 . 
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Figure 3: The lower bounds on m u + md as a function of s. These are obtained by saturating 
the spectral function with the pion pole contribution, and from the "ratio" method described 
in the text, s is in GeV 2 and m u + rrid in GeV. 
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Figure 4: Four examples of the hadronic spectral function, assuming different resonance 
structure and point of matching to the perturbative solution. Units are as in Fig. 2. The 
locations and widths of the resonances used are given in the text. The normalization A and 
the relative strengths q for the four cases are given in Table 1, along with the values of So 
and m u + m d used to derive the perturbative solution. 
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Figure 5: The hadronic spectral function assuming four resonances with the quantum num- 
bers of the pion. Units are as in Fig. 2. The locations and widths of the resonances are given 
in the text. The normalization A and relative strengths q for the three cases are given in 
Table 2. These have been adjusted to make the hadronic form join smoothly to the duality 
ansatz at s = 10 GeV 2 . The solid line corresponds to Case 1, the dotted line to Case 2, 
and the dashed line to Case 3. 
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Figure 6: The output values of m u + m& (run to fi = 1 GeV) for the three cases of the 
spectral function shown in Fig. 5. s is in GeV 2 and m u + in GeV. 
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Figure 7: Plots of the spectral functions p p qcd and phadronic used by JM/CPS [5, 4]. The scale 
of matching between the pQCD and hadronic solution is sq = 6.0 GeV 2 . To highlight the fact 
that the phadronic is dominated by the resonance contribution we also show pxn, i-e. Phadronic 
without the Breit-Wigner modulation factor. For convenience we plot p * 10 5 /(m s — m u ) 2 , 



so the units along the y axis are GeV 2 . 
been taken from Ref. [5]. 



The values, m s = 189 MeV and m u = 5 MeV have 
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